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Abstract 

In this paper, inspired by the fractional Brownian sheet, we first introduce the 
operator fractional Brownian sheet. Then, based on the martingale differences, we 
give a weak approximation of the operator fractional Brownian sheet. 

Keywords: Fractional Brownian motion, operator fractional Brownian motion, frac¬ 
tional Brownian sheet, operator fractional Brownian sheet, martingale differences. 

1. Introduction 

Self-similar processes, first studied rigorously by Lamperti m under the name ’’semi¬ 
stable”, are stochastic processes that are invariant in distribution under suitable scaling 
of time and space. There has been an extensive literature on self-similar processes. We 
refer to Vervaat [20] for general properties, to Samorodnitsky and Taqqu [Hj [Chaps.7 and 
8] for studies on Gaussian and stable self-similar processes and random fields. 

The fractional Brownian motion (fBm) as a well-known self-similar process has been 
studied extensively. Many results about weak approximation to fBms have been estab¬ 
lished recently. See [TUI [16] and the references therein. Due to some limitations of fBms, 
many researchers have studied the generalization of fBms. The fractional Brownian sheet 
introduced by Karnont m is one generalization of fBms. Scholars also have studied the 
limit theorems for this kind of processes. For more information, refer to mm and mm- 
The definition of self-similarity has been extended to allow scaling by linear operators 
on multidimensional space M. d , and the corresponding processes are called operator self¬ 
similar processes. We refer to m , m , m and the references therein. We note that 
Didier and Pipiras muni introduced the operator fractional Brownian motions (OFBM 
in short) as an extension of fBms and studied their properties. Similar to fBms, weak limit 
theorems for OFBMs have also attracted a lot of interest. Recently, Dai and his coauthors 
m-m presented some weak limit theorems for some kinds of OFBM processes. 

On the other hand, we should point out that, contrast to the extensive study on the 
extension of fBms, there is little work to study the extension of OFBMs. Inspired by the 
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study of the fractional Brownian sheet, we will introduce a new process, which we call the 
operator fractional Brownian sheet, and present a weak limit theorem for it. 

The rest of the paper is organized as follows. In Section 2, we introduce the operator 
fractional Brownian sheet and state the main result of this work. The proof of the main 
result will be detailed in Section 3. 

Most of the estimates of this paper contain unspecified constants. An unspecified posi¬ 
tive and finite constant will be denoted by C, which may not be the same in each occur¬ 
rence. Sometimes we shall emphasize the dependence of these constants upon parameters. 


2. Preliminaries and Main Results 


In this section, we first introduce the operator fractional Brownian sheet. Let (fl, J~, P) 
be a complete and separable probability space and {J'sj/, ( s , t ) £ [0, T] x [0, 5]} be a family 
of sub-<7-£elds of T such that T s ,t ^ JFg i f° r an y t < t,s < s, where S and T are positive 
numbers. Moreover, we denote by A t tS X(t',s') the increment of X over the rectangle 
((t, s), (t , s )], that is, 

A t>s X(t',s) = X(t,s) - X(t, s') - X(t,s) + X(t, s), 
where (t, s) < (t , s'). 

Let cr{A) be the collection of all eigenvalues of a linear operator A on M d . Let 


Aj 4 = min{ReA : A £ cr (A)} and = max{ReA : A £ cr(A)}. (2-1) 


Moreover, for any x £ M d , x T denotes the transpose of x. We first recall the operator 
fractional Brownian motion of Riemann-Liouville type introduced by Dai [81. Let D be 
a linear operator on W d with ^ < Xd,Ad < 1. For any t € M + , we define the operator 
fractional Brownian motion of Riemann-Liouville type X = {X(t)} by 

X(t)= [ ( t-u) D ~ I/2 dW(u ), (2.2) 

Jo 

where W{u) = { W l (u ),..., W d (u)} T is a standard d-dimensional Brownian motion. 

In order to introduce the operator fractional Brownian sheet, we need to review the 
relation between the fBm and the fractional Brownian sheet. By Alos et al. [T], the fBm 
B h of Hurst parameter H £ (0,1) has the following integral representation with respect 
to the standard Brownian motion B: 

B H {t)= f K H (t,a)dB(s), t> 0, 

Jo 

where Ku is the kernel defined on the set {0 < s < t} and it is given by 

K H {t,s) = C H {t- s) H ~^ + C H {\- H) [ (u - - (-)^ H )du, 

2 J s u 


with the normalizing constant Ch > 0 being given by 


Ch 


2oT(f -a) | 2 

T(i + a)T(2-2«) J 
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Recall that a fractional Brownian sheet with parameters a G (0,1) and (5 G (0,1) admits 
an integral representation of the following form (Bardina et al. [3]): 

W a ’P(t,s)=[ I K a {t,v)Kp{s,u)B(dv,du ), (t, s) G [0, T] x [0, S] (2.3) 
Jo Jo 


where B is a standard Brownian sheet and the deterministic kernel K a and Kp are the 
same as the Kh with H being replaced by a and /3, respectively. 

Inspired by the above, we can define the operator fractional Brownian sheet X = 
{X(t, s), (t, s) G [0, T] x [0, S']} as follows. 

Definition 2.1 The operator fractional Brownian sheet X = {X(t, s), (t, s) G [0, T] x 
[0, S]} is defined by 


X(t,s) = f f (t — u) 2 4( s — -y) 2 4 B(du,dv), (2.4) 

Jo Jo 

where B(du,dv) = (B 1 (du, dv ),..., B d (du, dv)) T with B l being independent copies of B, 
and | < A d, A d < 1- 

It is obvious that the equation (I2.4[) is well defined. Next, we study some properties of 
the process X. We first introduce the following notation. Let ||x ||2 denotes the usual 
Euclidean norm of x G M rf . For any linear operator A on R^, let ||t 4|| = maxiu^ ||Ax ||2 
denote the operator norm of A. Hence, we have 

Lemma 2.1 The random field X = {X(t, s), (t, s) G [0, T] x [0, S]} is an o.s.s Gaussian 
random field with exponent D. Moreover, X has a version with continuous sample path 
a.s.. 

Proof: We first check the operator self-similarity. For every c > 0, we have 


X(ct,cs ) = f [ (ct 
Jo Jo 

1 f ct f cs . u d_i v d_i 

2 / (t - ) 2 4 ( s - ) 2 4 dB(u,v) 

Jo Jo c c 


= c 




— u)+ 4 (cs — v)f 4 dB(u,v ) 

ret res 


= c D X(t, s), 


since 


B(cu, cv) = c^B(u, v ), (2.5) 

where X = Y means the processes X and Y have the same finite dimensional distributions. 

Next, we check the sample continuity. In fact, for any t, s G [0, T] x [0, S] with 
||t — s ||2 < 1, by some calculations, we have 


A s X(t) = 

fT fS / n_ — L R — L\ / R — L 

Jo Jo y {tl ~ u) + 3 


(s 2 — v)£ 4 ^jB(du, dv\2.6) 
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Hence, 


rT r S d 


A s x (t)\\l = Y] ( [ [ V ]F itj (t,s,u,v)B : >(du, 

i =i V/o 7o i=i 


dv 


(2.7) 


where 

F(i, s,u,u) = ((ti-n) P 4 -(si-v) 


D I D _ J 

2 4 


*2-u)+ 4 -(«2 -v)+ 4 ) = (Fi,j(t,S,U,v)) dxd . 


Noting /q T / 0 'Y^'j=\Fij(t, s , u , v )Bi{d u ,dv) is a Gaussian random filed, it follows from 
(12.71) that 


E 


A a X(t)||*l <C[ [ T [ S \\F(t,s,u,v)\ 
J L io Jo 


"dudv 


On the other hand, we have 

||F(t,s,u,u)|| < C\\Fi(t, s, u, v)\\ 2 x \\F 2 {t,s,u,v)\\ 2 , 

where 


and 


Fi(t,s,u,v) = (ti-u)* 4 -{s 1 -u) + 
F 2 (t,s,u,v) = (t 2 - u)_p 4 -(s 2 -v) 


D_ _ I_ 

2 4 


D_ _7 

2 4 


Now, we look at 

f T 2 

/ ||Fi(t, s, it, u)|| dudv. 

Jo 

Using the same method as in Dai, Hu and Lee [8], we have 

l \\{ti-u)l 4 -(si-u)l 4 \\ 2 du < C||(ti - si) T_? || 2 (U - Si). 
Jo 


Similar to (12.101) . 


[ \\F 2 (t,s,u,v)\\ 2 dv < C\\(t 2 - s 2 ) 2 4 || 2 (^2 — « 2 )- 

J 0 

From Maejima and Mason [18), (12.81) . (12.91) . (12.101) and (12.111) . we can get that 


E 


|A s Y(f)| 


< C 


(ti - s 1 ) Xd+ ^~ 25 x (t 2 - s 2 ) Ao+ 2 _2<S 


< C||i- 


i(A D + i-25)fc 


( 2 . 8 ) 

(2.9) 


( 2 . 10 ) 


( 2 . 11 ) 


( 2 . 12 ) 

The sample continuity follows from Garsia [6] and ()2. 121) . □ 

The aim of this paper is to prove a weak convergence to the operator fractional Brownian 
sheet via martingale differences. In order to reach it, we first recall some facts about 
martingale differences. Similar to Wang, Yan and Yu |2T], we will use the definitions and 
notations introduced in the basic work of Cairoli and Walsh [5] on stochastic calculus 
in the plane. For any 2 = (t,s) € [0, T] x [0, S ], let T z := Ft,T V J~S,s, the a— fields 
generated by F^t and J~s, s - Now, we recall the definition of the strong martingale. 
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Definition 2.2 An integrable process Y = {Y(z),z € [0, T] x [0, S]} is called a strong 
martingale if: 


(i) Y is adapted; 

(ii) Y vanishes on the axes; 

(iii) ¥.[^ z Y{z')\F z ] = 0 for any z < z' with the usual partial order. 

Let , F- j^)} ngN be a sequence such that for all 


E 


(n) 




F 


(n) 


= 0 , 


where F ^^ = F^ 


V j w bh F- 7 ^ and F^} being the a— fields generated by ^\ n } and 


respectively. Then we call ■ F ^)} nG pj a martingale differences sequence. 

It is well known that if the martingale differences sequence {^ n i} satisfies the following 
condition 


[nij [ns\ 


EE tel ?) 2 

2—1 2—1 

in the sense of L 1 , then the sequence 

\nt\ [nsj 

EE? 


t ■ s 


( n) 

,3 


2 — 1 2 — 1 


converges weakly to the Brownian sheet, as n goes to infinity (see for example, Morkvenas 
PCS].) Recently, Wang, Yan and Yu m extended this work to the fractional Brown¬ 
ian sheet. If {£!”'} is a square integrable martingale differences sequence stratifying the 
following two conditions 

lim n(£■”•*) = 1, a.s. (2.13) 

n—>-oo 

for any 1 < i, j < n, and 

max \^ n )\ <—,a.s. (2.14) 

1 <i,j<n n 


for some C > 1, then they [21] constructed a sequence { Z n } to converge weakly to the 
fractional Brownian sheet. Inspired by these work, we want to study the weak limit 
theorem for the operator fractional Brownian sheet we introduced in Definition 2.1. 
Define 


and 


JF) _ iAF) A n ) \T 

•i,j ***’> 2 ,j,d' ’ 


[nij 

= ij’ 

i=l j =1 


(2.15) 


(2.16) 
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where k = 1, • • • , d are independent copies of ■ 

From the above, we obtain that J r }^)} ngN is still a sequence of square integrable 

martingale differences on the probability space (il, J 7 , P). 

For any n > 1, (t, s ) £ [0,1] x [0,1], define 


X n {t,s) = 


/> 


— u)? 4 (s — v)? 4 B n (du,dv) 


D_ _ I_ 

2 4 


|_ntj |nsj 

= n2 E E 'ii'j 

i= 1 1=1 


(f — n) _? 4 (s — v) 


f i—1 / ,7 — 1 


D_ _ I_ 

1 4 dudv. 


(2.17) 


The main result of this paper is the following. 

Theorem 2.1 The sequence of processes {X n (t, s), (t, s) £ [0,1] x [0,1]} given by (12.171) . 
as n tends to infinity, converges weakly to the operator fractional Brownian sheet {X(t, s); (t, s ) £ 
[0, 1] x [0, 1]} given by (12.4|) . 

3. Proofs 

The proofs are based on a series of technical results. 


Lemma 3.1 Let {X n (t, s)} be the family of processes defined by (12.171) . Then for any 
s < t < u, there exists a constant C such that 


E 


|A s X n (t)|| 2 ||A t X n (u)|| 2 < C(u 2 - s 2 ) h {ui - si) H , 


(3.1) 


where H = Xp — 25 + \ with 5 > 0. 

Proof: We choose (t,s) < ( t',s ') £ [0, 1] x [0, 1]. Then 


X t , s X n (t ,s ) 

pt' rs' , D_l 

= Jt J V ~ U ^ 4 “ 

[ nt'j [ns'] 


D_ _ J 
2 4 


D_ _ I_ 

2 4 


s' — v)f 4 — (s — v)f 4 )B n (du,dv ) 

•- 1 — 1 — 

- v)f~^ - (^ - v)f ~ L4 )dudvrj$. 


i =1 1=1 
[rts'J 
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It follows from (12.141) that 


E 


\A t:S X n (t,s')\\ 2 2 


[nt' J |ns'J 

EE 

*=1 3 =1 


n 


n 



r i—i J j —i 

n n 


I V n n / 


[ns'J 

n 
. [nt' J 

<c(E" 

^ »=i 

|ns'J „i 

• E n 

J=1 


- «)+ 4 - ( 




n 


»-i. n n 


du 


124 _(JH 1 

j-i " ~ n n 


dv 


Using the Cauchy-Schwarz inequality, the above term can be bounded by 


c£,i‘i J jL ik 4^ - -j - (ifi - ~ i ii 2 *i 

n 

£}r,' J sL ii(¥ - »)f" i - (¥ - AAi 2 * 

n 

< C/l' ll( 4 r- - - (¥ - «)?“* II 2 * 


/o 


[ns'j 


“0 + 4||2 dw 


<cfo ik¥ -<•)-? 1 - (¥ -'ii 2 * 


/„ 1 ii(¥-*or i -(¥-nr , ii 2 *- 

From Dai, Hu and Lee [9], we obtain that 

r ii ^ ~ nt 

Jo 


D_ _ I_ 

2 4 — {s' — u) 


D I 

2 4 II2 


du < C(t' — s') H 


where H = Xp> — 25 + \- Then f|3.2|> can be bounded by 

r<(i n t'} ~ i n t] \ H f [ ns'] — [ns] 

C A n ) V n 

Hence, for any any s < t < u £ [0, 1] x [0, 1], we have 


H 


E 


||A s X(t)|| 2 ||A t X(n)|| 2 <C[E[A t X(u)||I]*[E[A,Jr(t)||I] 


2l| 


< c 


i} 


[nti] - [nsi] \ f / [nt 2 \ - [ns 2 ] \ f / [nui] - [nti] \ f / [nu 2 \ - [nt 2 ] \ f 


n 


n 


n 


n 


< c ,' [nu 1 \ - [nsi] \ H / [nu 2 \ - [ns 2 h # 


n 


n 


Hence, if n 2 — s 2 > then 


[nn 2 ] - [ns 2 ] 


n 


H 


< C\(u 2 - s 2 )\ h . 


(3.2) 


(3.3) 


(3.4) 


(3.5) 
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Conversely, if U 2 — S 2 < then either and £2 or t 2 and S 2 belong to a same subinterval 
[—, for some integer m. Hence (13.51) still holds. The second term follows a similar 

discussion. From the above argument, the proof is now completed. □ 

Using the criterion given by Bickel and Wickel [2], and noting that X n (t, s ) are null on 
the axes, and by lemma 3.1, we can get the following lemma. 

Lemma 3.2 The sequence {X n (t, s); (t, s) £ [0, 1] x [0, 1]} is tight. 

Now, in order to prove Theorem 12.11 it suffices to show the following theorem which 
states that the law of all possible weak limits is the law of a operator fractional Brownian 
sheet. 

Theorem 3.1 The family of processes X n (t,s ) defined by (12.171) converges, as n tends 
to infinity, to the operator fractional Brownian sheet X in the sense of finite-dimensional 
distribution. 

In order to prove Theorem 13.11 we need a technical result. Before we present this lemma, 
we first introduce the following notation. 

= (. kij(t,u)) dxd 


and 


[nt\ 


n 


D_ _ I_ 

2 4 




dxd' 


Then, we have 

Lemma 3.3 For any (t k , s k ), (t;, si) £ [0,1] x [0,1] and q,m £ { 1, ■ ■ ■ , d}, we have 
n 4 E"=i E"= 1 Ilk Ilk u)K^ q (s k , v)dudv 

n n 

fkl Ilk ^lm{^u)Kf nq {s h v)dudv{^k q ) 2 , 


n n 


converges to 


1 r 1 


(tk,u)K mtq (sk, v)K qtrn (ti,u)K m , q (si,v)dudv, a.s. 


'0 Jo 

as n tends to infinity. 

Proof: It is obvious that (13.61) is equivalent to 


E”=l n Ik 1 ^q,m( t k,u)du 1 K£ m (t h u)du 

n n 

• Ei=l n flk Kkq( S k, V ) dv ffk K£ q(si, V)dv(fff ) 2 . 


(3.6) 


(3.7) 


(3.8) 


Using the same method as the proof of Lemma 8 in Dai, Hu and Lee[9], we can prove the 
lemma. □ 

Next, we prove the Theorem 13.11 
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Proof of Theorem 13.11 Let a±, ..., oq G R and (t\, si), (tQ, sq ) € [0,1] x [0,1]. We 
see that the random variable 


Q 

Yn ^ ^ &k) 

k= 1 

converges in distribution, as n tends to infinity, to the Gaussian random vector 

Q 

X ^ ^ A: ■ ■Sfc)- 

fc=l 

By the well-known Cramer-wold device, for example see Whitt [23], in order to prove the 
above statement, we only need to show that as n —> oo 

bY n 4 bX, (3.9) 


where b = (b\, ■ ■ ■ , bg) and —> denotes convergence in distribution. 

For conciseness of the paper, let 

{t-u)l 4 (s-u)^ 4 = K(t,s,u,v) = {Ki(t,s,u, ■»),••• ,K d (t,s,u,v)} T , 


where 

Kj(t, s, u, v) = (Kj tl (t,s,u,v),- ■ ■ ,K jtd (t,s,u,v )). 

Then, we have 

d Q \nt\ |nsj 


»» = EEEE 


*_L « 2. 

n In 


n 


a k b q K q ( ^- ntk ^ , ,u.v)r$dudv 


, , ..-'q ± ^q\ 5 J 

/ j-i n n 


n n 


q= 1 /c=l 2—1 j = l 

d d Q \nt\ |rasj | , I I I 

E E E E E " 2 4 4 ^ 4 1 ' 

m=lg=lfc=lj=lj=l 


n n 


and 


d d Q r i „i 


“ = EEE 

m=l ^=1 /c=l 


0 ./O 


Q'kbqKq,m(j'ki ^4 V^B (dlL^dv). 


Since m = 1, ■ ■ ■ ,d are independent, in order to prove 113.91) . we only need to show 

d Q V nt \ L ns J 


EEEE 


r 2. 
n I n 


n ~ 4 i I j—l a k b Q K Q,rn(^^-> ^ u ^)^j,m dudv ^ 

q=lk=li=lj = l 

d Q ,1 „1 

EE/ / Q J kbqKq 1 m(j'ki v'jB (du,du). 

g=l fc=l ’ 


0 30 


(3.10) 
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For convenience, we introduce the following notation. 
d Q 

9=1 k= 1 

Then, ([3.1011 can be rewritten as 


( n ) 2 f n f 71 u ts ( \- n ^k\ L ns fcJ \A n ) J J 

- n / _! 
q= 1 k=l n n 


\nt\ |nsj d Q .i .i 

EE Y if ^^ 2^2 / a kb q K q , m (t k ,s k ,u,v)B m (du,dv). (3.11) 


*=i 1=1 


g=l fc=l 0 


Inspired by Wang, Yan and Yu [21], in order to prove (13.111) , we first prove the following 
Lindeberge condition 


n n 

„>“EE e K”Vi 

*=1 1=1 


{ly&W 


tt( n ) 

•'V1,1-1 


= 0 


(3.12) 


for all e > 0. 

In fact, we have 


J_ fk J_ 
n In 


n 



' i —1 J j —1 

n n 




n n 



[nt k ] [s4j 
n ’ n 

L?rtfcJ L s 4j 


n n 


, u, v)\dudv 
,u,v)\ 2 dudv. (3.13) 


We also note that for some 5 > 0 with A^> — i — 25 > 0, 


/ ||(f — it)? 4 || 2 du<C' / (1 — u)i° 2 2<5 du 

r i —1 / i —1 


< 


c hi- U )i D * 

do 


2 2< W 


(3.14) 


since (1 — u)i° 2 is decreasing in u. 

Noting the form of FT, it follows from (13.131) and (13.141) that 


,_L 

n In 


n- I I K q)m (2^1, ^, u ,v)g n ) m dudv) 
i—i /j-i n n ,J ' ‘ 


n n 


< Cn 2 ({i”! ra ) 2 «. 


(3.15) 


where 


5 n = [ (1 - it) 

Jo 


A_d —2<5 


du. 
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It follows from (13.131) and (13.151) that 


On the other hand, 


Y tt } ) 2 £ c E E 

q= 1 fc=l 


{|^ n) |> £ } = {|4 n) | 2 > £ 2 }. 


(3.16) 


Hence, from (13.171) and (13.161) . 

{ihl’l > e } £ {c’" 2 «Sh> 2<5 » > « 2 }- 

Consequently, 

E [(h; ) ) 2 i i|r <j )|>e) |41 J -i] < cE[ n 2 (4"> m ) 2 «;i {CnlKam)2JS>e2) |4r>,_ 1 ] 

- C ' 5 ” E [ 1 {Cn»Kg m )'«>c a ) 

for all i,j = 1,2,...,n. Hence, from (12.131) and (|3. 191) . 


(3.17) 


(3.18) 




(n) 


j-lj-U 


(3.19) 


EE®^) 2 * 




JO 


(n) 


i-lj-U 


*=1 j = 1 

n n 

~ EE C ' <5 n E [ 1 {Cn2(^ m )252 > eH 

2=1 j=l 


7 ; 


(n) 


i-lj-U 


— EE E [ 1 l c ' ,5 n>e 2 }] 0 (n -> oo), 

2=1 j = l 


because 5 n —> 0 implies 1 {Ci 52 > £ 2 }. —>• 0. 

In order to prove (13.91) . we also need to show that 

d Q r i r \ 


ioio i 2 ^ /*! 

EE[ y i n) ] ^ E [EE / / «fc^,2n(4,s fc ,u,u)H m (du,du)]. (3.20) 

L g=lfc=l^° "'° 


*=i i=i 

For convenience, we define 


H m (f, S, U, v) = E h q Kq,m{t, s, U, v). 
9=1 

Note that the right- hand of (|3.20l) is equivalent to 


12 /*! /*1 

i / B rn (ti, Si,u,v)B' n (tj,tj,u,v)dudv. 


2J=1 


0 J 0 


(3.21) 
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Next, we look at the left-side hand of (13.201) . In fact, we have 



(3.22) 


Hence 



k,l =1 


fim/NfeJ [ns k \ 

B (-,- ,u,v)dudv 

n n 


n 


n 



n n 


(3.23) 


n 


n 


Here, we also should point out that the entry K qjm (t,s,u,v) takes the form of 


d 


^ ' Kq,i (L ^)Kq i (S, (t, u)Ki m (s, v). 


Hence, it follows from Lemma 13.31 and (13.211) to (|3.23l) that (|3.20l) holds. 
From the above arguments, we can easily get that the theorem holds. 


□ 
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